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INTRODUCTION 


A  thick-walled  cylinder  is  used  for  a  variety  of  applications  in  the 
chemical,  nuclear,  and  armament  industries  where  large  internal  pressures  have 
to  be  withstood.  In  the  absence  of  residual  stresses,  the  cracks  usually  form 
at  the  bore  where  the  hoop  stress  developed  by  the  pressure  is  highest.  To 
prevent  such  failure  and  to  increase  the  pressure-carrying  capacity,  a  common 
practice  is  autofrettage  treatment  of  the  cylinder  prior  to  use  (ref  1).  This 
process  has  the  effect  of  producing  beneficial  compressive  residual  hoop 
stresses  near  the  bore  which  can  prolong  the  fatigue  life. 

The  determination  of  residual  stresses  in  autof rettaged  cylinders  has 
been  considered  by  many  investigators  using  different  mathematical  methods  and 
material  models  (refs  2-8).  Most  of  the  earlier  solutions  were  based  on  the 
assumption  that  the  material  behaves  elastically  on  the  release  of  the 
autofrettage  pressure.  However,  many  materials,  particularly  the  quenched  and 
tempered,  low  alloy  steels  generally  used  for  high  pressure  vessels,  exhibit  a 
significant  Bauschinger  effect  (ref  9).  In  a  recent  paper  (ref  10),  this 
author  presented  a  closed-form  solution  of  residual  stresses  in  autof rettaged 
tubes  based  on  a  theoretical  model  considering  the  Bauschinger  and  hardening 
effects  during  unloading,  but  neglecting  the  strain-hardening  effects  during 
loading.  A  more  general  theoretical  model  without  this  restriction  was  pro¬ 
posed  earlier,  but  only  part  of  the  final  results  were  shown  due  to  space 
limitation  (ref  11).  In  the  present  report,  the  complete  method  of  stress  and 
deformation  analysis  based  on  the  general  theoretical  model  is  stated.  The 
new  model  is  a  better  representation  of  the  actual  loading/unloading  behavior 

References  are  listed  at  the  end  of  this  report. 


in  a  high  strength  steel.  The  Bauschinger  effect  factor  is  treated  as  a  func¬ 
tion  of  overstrain.  The  strain-hardening  effect  is  taken  into  account  with 
different  parameters  used  for  loading  and  unloading  processes.  The  formulas 
for  calculating  stresses,  strains,  and  displacements  are  given  and  new  results 
of  residual  stresses  in  autof rettaged  high  pressure  vesssels  are  presented. 

MATERIAL  BEHAVIOR  AND  MODELING 

The  material  chosen  for  this  investigation  was  a  modified  4330  steel 
having  a  martensitic  structure.  A  description  of  its  chemical  composition  and 
various  heat  treatments  is  given  in  Reference  9  by  Milligan,  Koo,  and 
Davidson.  They  studied  behavior  by  utilizing  a  uniaxial  tension-compression 
specimen.  Figure  1  shows  the  stress-strain  curve  during  loading  and  unloading 
after  overstrains  in  tension.  The  stress-strain  curve  during  loading  can  be 
replaced  with  sufficient  accuracy  by  a  bilinear  elastic-plastic  model  shown  in 
the  figure  by  the  broken  lines.  For  the  plastic  portion,  the  yield  stress  (a) 
is  related  to  the  plastic  strain  ( c P )  by 

o/aQ  -  1  +  m</(l-m)  and  {  =  (E/cr0)cP  (1) 

where  Young's  modulus  (E),  tangent  modulus  (mE),  initial  yield  stress  (aQ) , 
and  the  Poisson's  ratio  (v)  are  the  material  constants.  The  nominal  yield 
strength  at  0.1  percent  offset  was  chosen  as  the  initial  yield  stress. 

Initially,  the  yield  stresses  in  tension  and  compression  are  approxi¬ 
mately  equal  so  that  the  material  can  be  considered  as  isotropic.  However, 
the  ratio  of  the  yield  stress  upon  reverse  yielding  to  the  initial  yield 
stress  is  strongly  affected  by  overstrain  in  a  high  strength  steel.  The 
values  of  the  Bauschinger  effect  factor  (BEF)  also  depend  on  the  0.1  percent 
offset  mentioned  above.  Taking  into  account  strain-hardening,  the  definition 


of  the  BEF  is 


BEF  =  (a0'-al)(/al  =  f(€lP) 


(2) 


where  aj,  ejP  are  the  yield  stress,  plastic  strain  just  before  unloading 
occurs,  and  a0'  is  tne  linear  drop  in  stress  until  reverse  yielding  begins. 
Figure  2  shows  the  Bauschinger  effect  factor  (f)  as  a  function  of  percent  ten¬ 
sile  overstrain  (cP).  The  graph  shows  a  decrease  of  the  BEF  with  increasing 
amount  of  tensile  prestrain  up  to  approximately  two  percent  at  which  point  it 
becomes  effectively  constant. 

In  most  of  the  plasticity  theories,  the  curve  of  reverse  loading  is 
uniquely  defined  by  the  curve  of  the  first  loading.  The  present  model  does 
not  assume  such  a  relationship.  The  experimental  stress-strain  curve  during 
unloading  is  used  directly.  A  piecewise  linear  representation  can  be  used, 
but  only  a  bilinear  approximation  was  chosen  here  as  shown  in  Figure  1.  In 
the  same  figure,  we  also  show  two  other  theoretical  models  in  dotted  lines  for 
this  material  subjected  to  unloading  with  further  plastic  flow.  They  are 
isotropic  and  kinematic  hardening  models.  A  summary  of  the  plasticity 
theories  given  by  Armen  (ref  12)  illustrates  the  fact  that  all  of  the  theories 
are  capable  of  treating  the  monotonic  loading  situation.  For  cyclic  loading, 
including  reverse  yielding,  kinematic  and  isotropic  hardening  represent  the 
limits  to  the  actual  behavior,  whereas  the  remaining  theories  are  falling 
anywhere  within  these  limits.  The  above  conclusion  indicates  that  none  of  the 
existing  models  can  accurately  represent  the  actual  material  behavior  as  shown 
in  Figure  1.  The  reverse  yielding  begins  at  a  stress  whose  magnitude  is 
smaller  than  that  predicted  by  all  existing  theories  and  the  slope  of  strain¬ 
hardening  after  reverse  yielding  is  much  larger  than  that  during  the  first 
loading  process.  The  new  bilinear  model  presented  here  is  a  better  approxima¬ 
tion  to  the  actual  material  behavior. 
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Figure  1.  Stress-strain  curve  during  loading  and  unloading. 

Choosing  a  new  coordinate  system  (a',e')  with  the  origin  at  the  point 
before  unloading,  we  have  for  the  plastic  portion  of  the  reverse  yielding 
curve 

o'/Oq  =  cr0'/a0  +  m’C'/d-m')  and  {'  =  (E/a0)c'P  (3) 

where  m'E  is  the  slope  of  the  reverse  yielding  curve  and  e'P  is  the  additional 
plastic  •‘train  during  unloading.  According  to  Eqs.  (1)  and  (2),  a0’  can  be 
expressed  as  a  function  of  plastic  strain  (ciP)  just  prior  to  unloading  by 

00‘/00  =  [1  +  md/d-mnu  +  f  ( Cl )  ]  *  g  ( Cl )  (4) 


With  experimental  data  f(eP)  shown  in  Figure  2,  we  can  calculate  a0'  and 
determine  the  initiation  of  reverse  yielding.  The  slope  of  reverse  yielding 
(m'E)  could  be  expressed  also  as  a  function  of  prestrain  (eP  and  e'P)  if  more 
experimental  unloading  curves  were  available.  Based  on  limited  information, 
the  slope  of  reverse  yielding  (m'E)  is  estimated  to  be  in  the  range  of  (9.4  to 
8.3)  x  10*  psi.  It  seems  that  this  slope  is  not  very  sensitive  to  prestrain 
and  a  constant  value,  say  m'  =  0.3,  may  be  assumed. 


Figure  2.  Bauschinger  effect  factor  as  a  function  of  overstrain. 

ELASTIC-PLASTIC  LOADING 

Consider  a  thick-walled  cylinder,  internal  radius  a  and  external  radius 
b,  which  is  subjected  to  internal  pressure  p.  The  material  is  assumed  to  be 
elastic-plastic,  obeying  the  Tresca's  yield  criterion,  the  associated  flow 
theory,  and  a  linear  strain-hardening  rule.  The  elastic-plastic  solution 
during  loading  has  been  found  by  Bland  (ref  3).  The  expressions  for  the 
radial  and  hoop  stresses  are.- 


(5) 


or/aQ  1  p2  l  p2  p 

«  -  (¥  1  +--)¥-  02  (~“  -  1)  ~  (l-fl2H°9  “  for  a  <  r  <  p 

Oq/o  o  2  b2  2  r2  r  (6) 

subject  to  oq  >  az  >  ar 

<*r/ao  1  p2  p2  (7 

=  -  (--  +  — )  for  p  <  r  $  b 

oe/aQ  2  b*  r*  (8 

ao  -  v(°r+ae)/ao  *  Eez/a o  (9 


E  u  ar  p2  E 

--  -  =  (l-2u) (l+v)  --  +  (l-v2)  --  -  v  --  ez  (10 

orQ  r  o0  r2  <r0 


and 

(£/°o)tz  =  {p/a°)  (n 

b2/a2-l 


where  p  =  0  (open-end),  1  (closed-end),  and  p  is  the  elastic-plastic  boundary 


relating  to  the  internal  pressure  p  by 


1  1 

p/o0  =  -  (l-p2/b2)  +  (l-02)1og(p/a)  +  -  02(p2/a2-l)  (12 


The  equivalent  plastic  strain  can  be  calculated  by 

(E/a0)eP  =  (  =  0! (p2/r2-l )  ,  in  (a  (  ip)  (13 

and 

1-m 

01  = - - -  ,  02  =  m0]/(l~m)  (14 

*3  (1-m) 

m  +  - - 

2  (l-v2) 


REVERSE  YIELOING 

If  the  pressure  p  given  by  Eq.  (12)  is  subsequently  removed  completely 
with  no  reverse  yielding,  the  unloading  is  entirely  elastic  and  the  solution 


is  given  by 


b*/a*-l 


r* 


(16) 


wrap  To  W  W  i  A  H  T  TO  A  m 


«r' 

ae’ 


a2’  =  u(ar*+a0.)  +  Eez'  (17) 

Eez’  =  -  (M-2v)p/(bVa*-l)  (18) 

Eu'/r  =  -[(l-v-nv)  +  (l+v>)b*/r*]p/(b*/a*-l)  (19) 

Let  a  double  prime  denote  a  component  in  the  residual  state,  i.e.,  a0"  = 
Oq  +  a0'.  Assuming  a  reduced  compressive  yield  strength  as  a  result  of  the 
Bauschinger  effect  and  using  Tresca's  yield  criterion  subject  to  or"  >  az"  > 
<r0",  the  reverse  yielding  will  not  occur  if 

<rr"  -  a q"  S  or"  =  fcr0[l  +  mC/( 1-m) ]  (20) 

Substituting  the  loading  and  unloading  solutions  into  Eq.  (20),  we  can  deter¬ 
mine  the  minimum  pressure  (pm)  for  reverse  yielding  to  occur.  The  equation 
for  pm  is  given  by 


Pm/^o  =  ~  (l-a«/b*)[l  ♦  ">fm/(l-"»)ni  +  f(tmH  (21) 


and  Cm  =  01 (Pm2/az-l ) . 

Equating  Eqs.  (12)  to  (21),  we  can  determine  the  minimum  amount  of 
overstrain  (pm)  and  calculate  the  minimum  pressure  (pm)  required  for  reverse 
yielding . 

ELASTIC-PLASTIC  UNLOADING 

Now  suppose  that  the  loading  has  been  such  that  the  internal  pressure  is 
larger  than  pm  given  by  Eq.  (21).  On  unloading,  yielding  will  occur  for  a  <  r 
*  p'  with  p'  >  p.  In  the  reverse  yielding  zone,  the  stresses  in  the  residual 
state  satisfy 

or"  -  o0"  =  a"  (22) 


7 


•  •  • 


assuming  that  ar“  >  az"  >  a0".  Since  the  yield  criterion  during  loading  is 

oq  -  or  =  o  (23) 

we  have 

<7r'  -  Oq'  ■  (Or“-Og")  +  [OQ-Or)  =  o"  +  O  =  O'  (24) 

where  a  and  o'  are  defined  by  Eqs.  (1)  and  (3),  respectively. 

The  material  is  assumed  to  obey  the  associated  flow  theory  and  a  linear 
strain-hardening  rule  during  unloading.  Following  the  procedure  in  Bland's 
work  (ref  3),  an  analytical  solution  for  elastic-plastic  unloading  can  be 
found.  The  stresses  in  the  reverse  yielding  zone  (a  (  r  (  p')  are  given 
by 

or'/aQ  =  P/o0  -  r(Cr.Ca)  (25) 

o0'/oo  =  or'/o0  -  g(0  -  m'C’/d-m’)  (26) 

C’  =  #i [g ( Cp ' )p' */r2  -  g ( C ) )  (27) 

where 

^  ( 1-m' ) 

/5l-  =  (l-m*)/[m’  +  --  >  02 '  =  m'P2'/(l-m')  (28) 

1 

HCr.Ca)  =  "  -  (02  '  /01 )  (P'/p) 2  (Cr"(a)9(Cp' )  ♦  (1-021 )G(Cr,Ca)  (29) 

fCa  1 

G(Cr-Ca)  =  Jcr  -  (f+fll>",g(C)df  (30) 

9(0  is  given  by  Eq.  (4)  and  Cr,  £a  can  be  calculated  by  Eq.  (12).  For  some 
special  cases,  e.g.  Reference  (10),  explicit  expressions  for  G(fr,(a)  can  be 
obtained.  In  general,  numerical  integration  is  needed. 

The  stresses  in  the  elastic  zone  (p'  <  r  4  b)  are 

1  (31) 
=  -  9 (Cp'HMp'/r)*  -  ( p ’ /b ) 2 ) 

2 


®r  * / 

'  /<*0 


(32 


The  location  of  the  reverse  yielding  zone  can  be  found  by  the  continuity  of 
or'  at  p'  from  Eqs.  (25)  and  (31). 

The  other  expressions  in  the  entire  tube  (a  <  r  <  b)  are 

°i'/o0  =  v(ar'+ae' )/oQ  +  Eez'/a0  (33) 

Bez'/a0  *  '  (M-2»)(P/o0)/(bVa*-l)  (34) 

(E/a0)u'/r  =  (l-2i>)  (1+v)  (ar' /aQ)  -  i>E ez’/o0 

-  (l-yMtP'/rJ^jCp*)  (35) 

The  residual  stresses  and  the  residual  displacement  are  found  by  addition 
ar"  =  ar  +  or‘  ,  o0"  =  oq  +  oq'  ,  oz"  =  oz  +  oz’ 

and 

u"  =  u  +  u'  (36) 

NUMERICAL  RESULTS  AND  DISCUSSION 

The  numerical  results  for  three  closed-end  thick-walled  cylinders  with 
various  wall  ratios  have  been  obtained.  The  material  constants  in  all  cases 

were  v  =  0.3,  E/o0  =  200,  and  m  =  0.01.  The  Bauschinger  effect  factor  (f)  was 

determined  during  computations  based  on  the  experimental  curve  f(eP)  as  shown 
in  Figure  8  of  Reference  9.  The  slope  of  unloading  after  reverse  yielding  was 
estimated  to  be  m'  =  0.3  for  a  high  strength  steel.  The  computations  for  the 
stresses,  strains,  and  displacements  during  elastic-plastic  loading  and  after 
complete  unloading  were  made,  and  some  of  the  numerical  results  are  presented 
graphically  in  Figures  3  through  10. 

Figure  3  shows  the  relation  between  the  pressure  factor  and  the  dimen¬ 
sionless  elastic-plastic  boundary  for  three  closed-end  tubes.  The  values  of 
the  pressure  factor  to  reach  100  percent  overstrain  are  0.543,  0.702,  0.891 
for  b/a  *  1.6,  2.0,  2.4,  respectively.  When  the  applied  pressure  in  a 
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partially  plastic  tube  is  removed,  residual  displacements  and  stresses  will 
occur.  Figure  4  shows  the  residual  displacements  at  the  bore  as  functions  of 
the  elastic-plastic  interface  for  wall  ratios  1.6,  2.0,  and  2.4.  By  combining 
the  graphical  results  shown  in  Figures  3  and  4,  we  can  obtain  the  functional 
relationship  between  the  applied  pressure  and  the  residual  displacement  at  the 
bore.  The  new  results  for  the  residual  displacement  shown  in  the  figure  are  a 
little  smaller  than  the  well-known  results  based  on  the  complete  elastic 
unloading  or  the  isotropic  hardening  material. 


Figure  3.  Internal  pressure  vs.  elastic-plastic  interface. 
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Figure  4.  Residual  displacement  vs.  elastic-plastic  interface. 

The  numerical  results  for  all  stress  components  during  loading  and  after 
unloading  have  been  obtained  for  tubes  with  various  wall  ratios.  Since  the 
radial  and  axial  stresses  are  not  as  important  as  the  hoop  stress,  only  the 
results  for  the  latter  are  presented.  Figure  5  shows  the  distribution  of  hoop 
stresses  before  and  after  unloading  in  a  closed-end  tube  of  wall  ratio  1.6 
corresponding  to  0  percent,  50  percent,  100  percent  overstrain,  i.e.,  p/a  = 
1.0,  1.3,  1.6,  respectively.  Reverse  yielding  will  not  occur  if  p/o0  <  0.46 
or  p/a  <  1.429.  For  a  100  percent  overstrain,  reverse  yielding  occurs  in  a 
small  region,  a  <  r  <  1.03a,  and  the  residual  hoop  stress  at  the  bore  is 
-0.474  cr0,  instead  of  -0.537  aQ  based  on  elastic  unloading.  Figure  6  shows 
the  residual  stress  distributions  in  an  autofrettaged  tube  of  wall  ratio  2.0 


for  p/a  =  1.2,  1.6,  and  2.0.  Reverse  yielding  will  not  occur  if  p/aQ  <  0.583 
or  p/a  =  1.373.  When  the  pressure  to  reach  100  percent  overstrain  has  been 
removed  completely,  reverse  yielding  occurs  in  a  C  r  <  1.173a  and  the  residual 
hoop  stress  at  the  bore  is  -0.377  aQ.  Additional  results  of  residual  hoop 
stresses  are  shown  in  Figure  7  for  an  autof rettaged  tube  of  wall  ratio  2.4 
with  p/a  =  1.2,  1.6,  2.0,  and  2.4.  Reverse  yielding  will  not  occur  if 
p/o0  <  0.647  or  p/a  <  1.358.  After  complete  unloading  from  100  percent 
overstrain,  reverse  yielding  occurs  in  a  (  r  ^  1.315a  and  the  residual  hoop 
stress  at  the  bore  is  -0.389  aQ.  As  can  be  seen  in  Figures  6  and  7,  the 
reverse  yielding  zone  becomes  larger  as  the  amount  of  overstrain  increases. 


Figure  5.  Stress  distribution  before  and  after  unloading  in  a  tube 
with  b/a  =  1.6. 
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Figure  7.  Residual  stress  distribution  in  an  autofrettaged  tube 
with  b/a  =  2.4. 

It  is  of  interest  to  compare  the  new  results  for  residual  hoop  stresses 
with  those  based  on  the  other  theoretical  models.  Figure  8  shows  a  comparison 
of  residual  hoop  stresses  based  on  three  theoretical  models  for  an  auto¬ 
frettaged  tube  of  b/a  =  p/a  =  2.0.  The  dotted  line  represents  the  well-known 
solution  based  on  the  assumption  of  complete  elastic  unloading.  The  same 
solution  is  obtained  when  either  the  isotropic  or  the  kinematic  hardening 
model  is  used.  The  broken  line  represents  the  solution  based  on  a  model 
including  the  Bauschinger  effect  but  with  the  same  hardening  parameter  m1  *  m 
*  0.01.  This  solution  is  close  to  Parker's  solution  neglecting  hardening  com- 
petely  with  m'  =  m  =  0  (ref  8).  As  can  be  seen  from  this  figure,  the  new 
result  represented  by  the  solid  line  is  quite  different  from  the  other  two 
solutions.  This  comparison  also  demonstrates  that  the  Bauschinger  effect  and 


strain-hardening  during  loading  and  unloading  can  have  a  significant  effect  on 
the  residual  stresses,  especially  near  the  bore. 


Figure  8.  A  comparison  of  three  theoretical  predictions  (b  =  p  =  2a). 

For  a  closed-end  tube  of  wall  ratio  2.4  subjected  to  100  percent 
overstrain,  we  show  the  comparison  of  residual  hoop  stresses  based  on  three 
models  in  Figure  9.  The  dotted  lines  represent  the  solutions  based  on  the 
isotropic  hardening  model  (ref  3).  Reverse  yielding  occurs  even  if  the 
Bauschinger  effect  is  neglected.  The  solid  and  broken  lines  represent  the  new 
model  and  the  model  with  m'  =  m,  respectively.  As  can  be  seen  from  these 
figures,  the  numerical  results  of  the  residual  hoop  stresses  based  on  three 


models  are  quite  different  in  the  reverse  yielding  zone.  The  new  result  is 
considered  more  accurate  for  our  problem  because  it  is  based  on  a  more  general 
theoretical  model  considering  the  Bauschinger  and  hardening  effects  during 
loading  and  unloading. 
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Figure  9.  A  comparison  of  three  theoretical  predictions  (b  =  p  =  2.4a). 

In  order  to  check  the  new  theoretical  prediction  of  residual  stresses,  a 
search  of  available  experimental  results  has  been  made.  Figure  10  shows  a 
comparison  with  the  results  based  on  the  x-ray  method  (ref  13)  and  the  ultra¬ 
sonic  method  (ref  14)  for  a  155  mm  autof rettaged  gun  tube.  A  good  agreement 
has  been  reached  in  this  comparison.  The  information  on  the  residual  stress 
distribution  has  been  used  in  calculating  the  stress  intensity  factors  (ref 


15)  and  predicting  the  fatigue  life  of  gun  barrels  (ref  16).  These  results 
will  be  published  in  the  near  future. 


(r-a )  /(b-a) 


Figure  10.  A  comparison  of  theoretical  prediction  and  experimental 
measurements  in  a  155  mm  autof rettaged  gun  tube. 

CONCLUSIONS 

A  new  theoretical  model  for  the  high  strength  gun  steel  has  been  pro¬ 
posed.  The  real  Bauschinger  effect  factor  has  been  used  to  determine  the 
range  of  elastic  unloading.  Different  rates  of  strain-hardening  during 
loading  and  unloading  have  been  taken  into  account. 

A  new  analytical  solution  for  calculating  the  residual  stresses  and 
strains  with  reversed  yielding  has  been  obtained.  The  numerical  results  for 
the  residual  stress  distributions  in  three  autof rettaged  gun  tubes  have  been 


given.  The  new  results  indicate  that  the  influence  of  the  Bauschinger  and 
hardening  effects  on  the  residual  stresses  is  significant. 

A  good  agreement  with  two  experimental  results  for  a  155  mm  autofrettaged 
gun  tube  has  been  reached.  The  new  method  is  reliable,  efficient,  and  econom¬ 
ical.  The  new  results  have  great  implication  to  the  fatigue  life  estimation 
of  gun  barrels. 
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